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We study the exchange constants of MnV^CU using magnetic force theorem and local spin density 
approximation of density functional theory supplemented with a correction due to on-site Hubbard 
interaction U. We obtain the exchanges for three different orbital orderings of the Vanadium atoms 
of the spinel. We then map the exchange constants to a Heisenberg model with single-ion anisotropy 
and solve for the spin-wave excitations in the non-collinear, low temperature phase of the spinel. 
The single-ion anisotropy parameters are obtained from an atomic multiplet exact-diagonalization 
program, taking into effect the crystal-field splitting and the spin-orbit coupling. We find good 
agreement between the spin waves of one of our orbital ordered setups with previously reported 
experimental spin waves as determined by neutron scattering. We can therefore determine the 
correct orbital order from various proposals that exist in the literature. 
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I. INTRODUCTION 

Transition metal oxides (TMO) are a class of solid- 
state materials that exhibit a rich variety of physi- 
cal phenomena 1 . Among them, magnetic cubic spinels 
AV2O4 have recently attracted much attention due to 
geometrically frustrated corner sharing tetrahedral net- 
work formed by the V atoms (also known as a py- 
rochlore lattice)^. An interesting example is represented 
by MnV2C>4 which is the spinel having additional mag- 
netic Mn ions. It exibits an orbital ordering (00) that 
occurs at finite T as a thermal phase transition: At room 
temperature, single crystalline MnV2 04 is a cubic para- 
magnet (PM) where Mn sites occupy the centers of oxy- 
gen tetrahedra (Mn04 units), while V sites occupy the 
centers of oxygen octahedra (V06 units) . As T is lowered 
there is a magnetic phase transition at Tp — 56 K from 
the high-T PM phase to a cubic ferrimagnetic (FEM) 
phase, with the Mn and V moments anti-aligned. When 
T is lowered further below T$ = 53 K there is a second 
transition to a tetragonal, non-collinear FEM with trig- 
onal distortion of V06 octahedra and orbital ordering of 
V 3+ 3d 2 electrons^. The structural transition and or- 
bital ordering occur together via the Jahn-Teller effect, 
which states that when there is an orbital degeneracy, 
as in the V atoms, structural distortion occurs to re- 
move the degeneracy and lower the total energy. The or- 
bital ordering is also accompanied by a reduction of the 

V magnetic moments due to the formation of the elec- 
tron orbital moment (finite orbital angular momentum). 
The orbital moment, m ~ 0.34, is anti-aligned with the 
spin moment, m s ss 1.65, giving the total moment of 
to rj 1.31®. The reduced value of V moment has been 
reproduced by an earlier first-principles work in Ref. HI 
and is explained by the spin-orbit coupling (SOC) on the 

V 3c? 2 which generally favors anti-alignment of spin and 
orbital angular momenta for T below the energy scale of 

socP. 

The local tetrahedral and octahedral coordination of 



the Mn and V sites, respectively, results in lifting the de- 
generacy of their five 3d orbitals known as crystal-field 
(CF) splitting. From group-theoretical consideration^, 
the 5-fold degenerate 3d level is split into a 2-fold de- 
generate e g and a 3-fold degenerate t 2g manifold. For 
tctrahedrally coordinated Mn, the e g is lower in energy 
as compared to the t 2g , whereas for octahedrally coor- 
dinated V, the e g — t 2g splitting is opposite. There is 
also Hund's rule splitting of up and down spins, which 
is greater than the CF splitting. In the stoichiometric 
crystalline environment, Mn has an outer shell high-spin 
5 = 5/2 configuration of 3d 5 and a valence of +2: all 
the up-spin 3d orbitals are occupied, and the down spin 
ones are empty. V has a valence of +3, an outer shell 
configuration of 3d 2 , and total spin S = 1: in this case, 
2 electrons must occupy the 3 t 2g xy, yz, zx orbitals. 
In the high temperature cubic phase, these latter three 
are degenerate, while in the low temperature tetragonal 
phase, where the unit cell is slightly compressed along the 
c-axis, the xy is lowered in energy while yz, zx remain 
degenerate. Thus, in the tetragonal (low-T) phase, one 
electron on V occupies the xy, and the second electron 
has the freedom to occupy either yz, or zx. The total 
angular momentum of the Mn atom is quenched, L = 0, 
due to the equal occupation of all the z-axis projections 
of L, which makes the total projection L z — 0. On the 
other hand, the total angular momentum of V is not fully 
quenched: The partial occupation of the yz and zx gives 
an effective total angular momentum L = 1 for V. The 
fact that L ^ implies that there maybe non-negligible 
effects of SOC in the V atoms 2 . Further, this is a hint 
that the yz and zx form complex linear combinations of 
one-electron states, since only such a complex state can 
carry non-zero angular momentum. The freedom of the 
second electron of V to occupy yz, zx, or some linear 
combination of the two gives rise to the possibility of 
long-range orbital order in the low temperature phase. 

The number of possibilities for the occupation of the 
second electron is infinite, due to the infinite different 
ways in which one can form linear combinations of the yz, 
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and zx orbitals. However, two simple choices exist and 
both have been studied theoretically in mean field mod- 
els. One is the Antiferro-Orbital Order (AFOO) with 
alternate occupation of the yz and zx along the c-axis, 
i.e.: the same orbital is occupied in a given a6-plane but 
the other orbit al is occupied in the adjacent planes above 
and below^ia as shown in Fig. [ija) . This order has the 
point -group symmetry IA\/a. The second is the Ferro- 
Orbital Order (FOO) where the same orbital is occupied 
on all V-atoms^, giving the symmetry group IAi/amd, 
as shown in Fig. [ljb). In particular, in the latter orbital 
order the linear combination can be complex, in which 
case there will be a non-zero orbital angular momentum 
and a magnetic moment associated with itP. Spin-orbit 
coupling can stabilize the finite orbital moment, since the 
energy is lower for anti-parallel alignment of L and S. 




FIG. 1. (a) Schematic illustration of the initial real antiferro- 
orbital order of the type I (AFOO I) with L = on the four 
corners of the V tetrahedron. The lower and upper horizontal 
bonds are in the ab plane. The red spheres are the V atoms. 
The lower ab plane has yz orbitals occupied on all V, while 
the upper ab plane has zx occupied on all V. (b) Schematic 
illustration of initial ferro orbital order on all four corners of 
the V tetrahedron where an electron occupies the same real 
linear combination of yz and zx on all V sites. Note that the 
self-consistent solution breaks this symmetry and results in 
an electron occupying alternately ip+ = (ip yz + 4>zx)/V2 and 
ip- — (V>j,z — i>zx)/ V% along the c-axis. We refer to this order 
as antiferro-orbital order of type II (AFOO II). The indices 
i, j, k, 1 denote the inequivalent V sites in the FCC primitive 
cell. 

In both of the above proposals, trigonal distortion of 
the V06 octahedra is not taken into account, but it is 
known to be large in MnV204 as compared to other 
vanadates. This distortion lifts the residual degeneracy 



between yz, zx, and the degeneracy of the e g manifold. 
In this case, the above 00 proposals are not necessarily 
correct as these assume degeneracy between the yz, zx 
orbitals. The trigonal d istort ion could possibly be due 
to the Jahn-Teller effectp2EQ anc j nas indeed been ob- 
served in a previous first-principles worlPthat used local 
spin density approximation (LSDA) of density functional 
theory (DFT^P^l supplemented by the correction due to 
on-site Hubbard interaction Ci^for correlation strengths 
U > 2 eV. In that work, in addition to a tetragonal re- 
laxation (compression) along the c-axis, structural relax- 
ation of the O positions is performed and a trigonal dis- 
tortion of the VOg octahedron with a concomitant low- 
ering of symmetry from IAi/amd to 74i/a is found. In 
this relaxed structure, due to the low symmetry at the 
V sites, all five 3d orbitals get mixed and the resulting 
five one-electron eigenstates of the on-site V block of the 
full LSDA+U Hamiltonian have different eigenenergies. 
By projecting the converged density onto an atomic or- 
bital basis using so called non-linear muffin-tin orbital 
(NMTO) downfoldingSl, the authors of Ref.Hfind a dif- 
ferent electron occupation order from the ones proposed 
above, namely, the first electron occupies the lowest en- 
ergy eigenstate, and the second occupies the next higher 
energy eigenstate, which is the same on all V sites, but 
rotated alternatively by 45° along the a6-chains due to 
the staggered trigonal distortion. Thus, the same orbitals 
are occupied on all V sites, akin to the ferro-orbital or- 
der, but nevertheless the symmetry is the /4i/a expected 
of antiferro-orbital order due to the trigonal distortion. 

The low-T magnetic excitations of the compound have 
been mapped along high -sym metry directions using in- 
elastic neutron scatterin g ^ 3 * 15 *. At the V point, these 
excitations are gapped for the acoustic modes, indicat- 
ing the presence of single-ion anisotropy, which essen- 
tially occurs due to the interplay between SOC and 
crystal-fields^. In Ref. [HI the authors start with a 
nearest-neighbor Heisenberg Hamiltonian including the 
anisotropy term and calculate spin-wave spectra and 
corresponding eigenmodes using linear spin-wave theory 
(LSWT) for the non-collinear, tetragonal phase. By fit- 
ting the spectrum to inelastic neutron scattering data 
(performed at sample temperature of 5 K), they were 
able to determine the exchange couplings between Mn- 
V, V-V in a&-plane, and V-V between a&-planes along 
the c-axis. They find all exchanges to be AFM with the 
following values: 

• Jmb-V = —2.82 meV 

• J v-v = - 9 - 89 meV 

• jc_ v = -3.08 meV 

The interplanar coupling between V atoms, Jy_ v , 
along the c-axis is unusually large for AFOO because 
such an alternate orbital occupation in the vertical di- 
rection would yield negligible orbital overlap, and would 
also be ferromagnetic (wrong sign) by the Goodenough- 
Kanamori rules^. The alternate proposal, FOO, would 
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be consistent with these results, but would have the 
wrong symmetry, IAi/amd. The symmetry group of 
this spinel vanadate has been established conclusively 
as I /a by a synchrotron x-ray study^ which supports 
AFOO, but contradicts with the large value of Jy-V- A 
possible resolution of this puzzle is that trigonal distor- 
tion has been ignored in these simple proposals. With 
trigonal distortion, we expect a more complex orbital 
ordering which has the requisite symmetry I4\/a and 
would give the observed (or fitted) J V -v a l° n g the c- 
axiiP^. This is exactly what has been found in the ab- 
initio work of Ref. H| Their physical picture has also 
received some support by a recent 51 V NMR work of 
Ref. [18] and by analytical model of Ref. [T9l 

In this work, we report our study of MnV204 based 
on the LSDA+U method and using linear muffin-tin or- 
bital (LMTO) basis set to solve the electronic structure 
problerrl 20 * 21 *. We calculate the pair-wise interatomic 
magnetic exchange interactions (J) between all mag- 
netic atoms using linear response theory and magnetic 
force theorerr P 2 * 23 ! , including the single-ion anisotropics 
(D) for Mn and V found by the exact diagonalization 
procedure^. We then use the obtained J and D as pa- 
rameters in a Heisenberg Hamiltonian with anisotropy 
to derive the spin-wave spectra in a semiclassical ap- 
proximation. We explore three initial orbital ordering 
scenarios: [1] Antiferro, [2] Ferro, and [3] Complex ferro 
+ SOC in the density matrix of the 3d shell of V to see 
how they affect the obtained exchange interactions. We 
also performed non-collinear magnetic electronic struc- 
ture calculations. Although we do not relax the struc- 
ture, and hence do not incorporate trigonal distortion, 
we find a good fit between our theoretical spin-waves for 
the SOC complex ferro-orbital order and the experimen- 
tal spin waves. We conclude that trigonal distortion is a 
higher-order consequence of the orbital order and of elec- 
tron correlations, and as such has a negligible effect on 
the exchange constants. Since LSDA+U has many solu- 
tions, we make the argument that our converged charge 
density maybe as good as the one found in Ref. |4] for the 
purposes of calculating the exchange interactions. We 
can further justify this point of view using the approach 
of Kugel-Khomskii (KK) models, in which both spin and 
orbital degrees of freedom are dynamical variables cou- 
pled to each other (these models are derived from multi- 
orbital Hubbard models using second-order perturbation 
theory)^. Simultaneous spin and orbital order is ob- 
tained by minimizing the mean-field total energy of the 
KK Hamiltonian with respect to the local orbital and spin 
order parameters. There is not yet any lattice distortion 
(indeed, the lattice is not a dynamic variable in the KK 
Hamiltonian), however, with a non-zero electron-phonon 
interaction (which is generically present in all materials) 
any orbital order (or polarization) will autom atically lead 
to a lattice distortion as a cons equenc S 10 ^ 24 ^. 

Along these lines, an example of orbitally driven lat- 
tice distortions is discussed in the context of the spinels 
MgTi2C>4 and CuIr2S4 by Khomskii and Mizokawa in 



Ref. [25] These authors consider orbital ordered Peierls 
distortion and spin singlet formation along chains of Ti 
and Ir atoms. The lattice distortion then follows as a 
consequence. 

Our paper is organized as follows. We begin with a 
discussion of the proposed orbital orders and their elec- 
tronic structures in Section II. We present our results for 
exchange interactions and comparisons with experiment 
in Section III. We end with the conclusions in Section IV. 



II. PROPOSED ORBITAL ORDERS AND 
THEIR ELECTRONIC STRUCTURES 

We have done LSDA+U calculations to model the elec- 
tronic structure for all three thermodynamic phases of 
MnV204. We describe our results in the following sub- 
sections for the T — phase only since this is the phase 
which exhibits orbital ordering and non-collinear mag- 
netism. Our results for the other finite-T phases may 
be found in Ref. US] For the magnetic phases we use 
the same values of U and J for both Mn and V corre- 
lated 3d shells. The use of the same U on Mn and V is 
justified because these elements have atomic numbers 25 
and 23, and are thus expected to have similar interac- 
tion strengths^. The Coulomb and exchange parameters 
in the solid state are generally screened, and hence re- 
duced by a considerable amount from their bare atomic 
values^. The structural parameters for all three phases 
are taken from experiment In the cubic phase, the lat- 
tice constant is 16.0746 a.u., and in the tetragonal phase 
it is 16.12 a.u. with a small tetragonal distortion ratio of 
£ = 0.98. 

a 

The non-collinear orbital ordered phase occurs when 
the temperature is reduced below Tg — 53 K. This phase 
transition results simultaneously in: [1] a structural tran- 
sition from cubic to tetragonal; [2] the canting of V mo- 
ments from a collinear ferrimagnetic (FEM) to a q = 
non-collinear FEM spin order with non-zero components 
in the a6-plane; and [3] a long-range orbital order in the 
V t2 g manifold. We model the electronic structure of 
this phase using LSDA+U method with [7 = 5 eV and 
J = 1 eV, but starting the self-consistency loop after 
imposition of the initial orbital order(s) in the Hubbard- 
U density matrix (further described below), along with 
tetragonal distortion and two different magnetic config- 
urations: [1] collinear, as in the intermediate phase, and 
[2] non-collinear, which is in fact the correct magnetic or- 
der for this phase. The converged charge density for the 
low-T collinear calculation was used as the initial charge 
density for the correct low-T non-collinear calculation. 
The orbital order that is finally obtained after reaching 
the self-consistency is taken to be the correct metastable 
solution within this approximation and specified initial 
condition(s). 

We initialize the V 3d density matrix to a particu- 
lar orbital order by specifying orbital occupation num- 
bers in the atomic basis. This means we initially spec- 
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ify only the diagonal components (occupation numbers) 
(n X yt)i ( n yz-f)> ( n zxt) °f the density matrix for all four 
V atoms' 3d shells and set the off-diagonal elements to 
zero. The full complex density matrix in the atomic ba- 
sis is (nmer,m'cr') (where m, m! and a, a' are the 3d orbital 
and spin indices respectively) and includes off-diagonal 
components as well. As a result of the electron-electron 
interactions, during the self-consistent cycle non-zero 
off-diagonal components of the density matrix develop 
(since the interactions mix the single-particle 3d orbitals 
at the Hartree-Fock (HF) mean-field level). This means 
the true occupied orbitals are some linear combination of 
the atomic basis functions. After convergence is reached, 
the final density matrix, which is no longer diagonal in 
the (mcr, m'cr') basis is diagonalized. The resulting eigen- 
vectors and eigenvalues give the 'correct' single-particle 
HF wave functions and their occupation numbers respec- 
tively. In the basis of these eigenfunctions, the density 
matrix is once again diagonal, and its non-zero entries 
signify the true orbitals which are occupied within the 
mean-field approximation of LSDA+U. We are thus able 
to identify the orbital ordering that results after conver- 
gence is attained. We describe the final orbital orders 
below for the collinear magnetic solutions only, since we 
find that the electronic structures, and therefore the or- 
bital orders, of the non-collinear configurations are not 
significantly different from the corresponding collinear 
ones as discussed below. 



A. Anti-Ferro Orbital Order I: IA\/a symmetry 

In the low-T orbitally ordered phase, the tetragonal 
distortion occurs to break the degeneracy of the t 2g in 
both V and Mn. There is no orbital freedom to place 
the electrons in the Mn 3d. In the V, the energy of 
xy gets lowered, so the first electron occupies xy. The 
second electron then has the freedom to occupy the re- 
maining degenerate orbitals yz or zx. Figure [T|a) shows 
the initial orbital occupations with /4i/a symmetry. In 
this scenario, the second electron of V occupies either yz 
or zx alternately along the c-axis (antiferro 00), and 
the same orbital within each ab plane (Each V chain 
within an ab plane has the same orbital occupied.) The fi- 
nal converged density matrix of the V 3d subspace shows 
that this initial order persists until convergence, although 
some small off-diagonal components do emerge. We label 
this order 'AFOO F. 

The collinear spin fat bands of V t 2g electrons are 
shown in Fig. [2]ja), and the same for non-collinear spins 
in Fig. [2jb). The occupations of the t 2g -\ bands, as 
shown by the partial characters, reflects the converged 
orbital order, as well as the FEM spin configuration. We 
also find that imposition of orbital order opens a gap 
of about i? gap = 1.67 eV at the Fermi level leading to 
an insulator state. The qualitative features of the band 
structure and partial characters does not change upon 
canting the V moments to the non-collinear configura- 
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FIG. 2. (a) The V t 2g -t bands for U = 5 eV in the low-T 
tetragonal phase with a collinear ferrimagnetic spin configu- 
ration and real antiferro-orbital order of the type I (AFOO 
I) as discussed in text, (b) Bands for the same setup as in 
(a) but with a non-collinear ferrimagnetic spin configuration. 
In both panels, the partial characters of the xy-^f, yz-f, zx-^[ 
orbitals are for the sublattice i, j V atoms. We see that due 
to the orbital ordering, the zx-f is occupied while the yz-~\ 
is mostly unoccupied. The occupations of these two orbitals 
are reversed for the sublattice k, 1 V atoms on the adjacent 
parallel ab planes along the c-axis. The sublattice indices are 
defined in Figs. [TJ[5] There is a band-gap of -E gap = 1-67 eV. 



tion: The band gap remains robust and the phase is still 
insulating. 



B. Anti-Ferro Orbital Order II: 74i/a symmetry 

The next simplest initial order has the second t 2g elec- 
tron occupying the same real linear combination of yz 
and zx on all V sites, with equal weight for both or- 
bitals, see Fig. [ljb). This initial order has I4i/amd 
symmetry. The real linear combination implies that the 
orbital angular momentum is zero, L = 0. We implement 
this by setting the initial mean occupations: (n xy ) = 1, 
( n yz) = (n Z x) = 1/2 (ferro 00), and the off-diagonal ele- 
ments to be zero. For this setup, the initial order does not 
persist until convergence is reached. Instead, there are 
significant non-zero off-diagonal elements, on the same 
order as the occupied diagonal elements, in the final den- 
sity matrix. Upon diagonalizing this final matrix, the 
orbital order we get has the second electron occupying al- 
ternately ip + = (ipy Z +ip zx )/V2a,ndip_ = {ip yz -ip zx )/V2 
along the c-axis, which again has the same IA\/a sym- 
metry considered in the preceding subsection. Thus, we 
start with an orbital order with I4i/amd symmetry, but 
the self-consistent solution breaks certain discrete sym- 
metries and results in an order with 74 1 /a symmetry. We 
thus label this order 'AFOO IF. 

The collinear spin fat bands of V t 2g and e g are shown 
in Fig. p3la), and the same for non-collinear spins in 
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FIG. 3. (a) The V t 2g ~t bands for U = 5 eV in the low-T 
tetragonal phase with a collinear ferrimagnetic spin configu- 
ration and real antiferro-orbital order of the type II (AFOO 
II) as discussed in text, (b) Bands for the same setup as in 
(a) but with a non- collinear ferrimagnetic spin configuration. 
In both panels, the partial characters of the xy-^f, yz-f, zx-^ 
orbitals are for the sublattice i, j V atoms. Note that the 
yz and zx partial characters have nearly identical dispersions 
due to their equal weight in the occupied orbital. There is 
again a band-gap for this orbital-order too. 

Fig. j3^b) . Qualitative features of the band structures do 
not change significantly between the collinear and non- 
collinear spin configurations. In both plots, the occupa- 
tions and dispersions of the yz and zx bands are nearly 
identical since these orbitals contribute equal weights to 
the true orbitals, although their relative signs in the lin- 
ear combinations might differ in these depending on the 
particular V atom. We also find an insulating band-gap, 
which in this case is smaller than for 'AFOO I.' 



C. Complex Ferro Orbital Order: Iii/amd 
symmetry 

The last 00 has one electron in xy as before, and 
the second electron in the spherical harmonic L z = — 1, 
S z = +1/2 state, which is a complex linear combination 
of yz and zx, on all V sites. This is an initial ferro-orbital 
order, but with SOC switched on and non-zero orbital 
angular momentum. The initial density matrix configu- 
ration persists until convergence. This calculation is car- 
ried out using LSDA+U+SO. This scenario is also illus- 
trated by Fig. [ljb), except that each V atom now carries 
a non zero orbital angular momentum of magnitude one 
due to the complex linear combination; hence, there is a 
uniform orbital order on all V atoms with L = 1 in the 
3d density matrix. The reason for choosing the opposite 
z-projections for L and S is that spin-orbit interaction 
lowers the energy for such a setup, compared to the case 
of having the same sign for both z-projections. 

In Fig.|4^a) we present the band structure of MnV204, 



FIG. 4. (a) The V t 2g -t band characters for the low-T tetrag- 
onal phase with a collinear FEM spin configuration and com- 
plex ferro-orbital order with spin-orbit coupling (SOC-FOO) 
as discussed in text, (b) Bands for the same setup as in (a) 
but with a non-collinear ferrimagnetic spin configuration. In 
both panels, the partial characters of the xy-f, yz-^[, yz-^ or- 
bitals are for the sublattice i, j, k, IV atoms. There is a 
band-gap of -E gap = 1.76 eV. The V xy-\, band lies above Ef, 
the zx-f bands coincide with the yz-f so we omit it, and fi- 
nally the zx-\. bands are in the same energy region as the yz-\. 
so again we omit it. 



for the collinear magnetic configuration, with the V t 2g -\ 
partial characters of the SOC uniform orbital order. We 
find that the V t^g-i and e g characters are above the Ep 
as expected. For Mn, all the 3d- 4- are below Ep, while 
the 3d-f are above. There is a band gap of 1.76 eV. 
Within LSDA+U a half metallic solution was found in 
Ref. HI with only the f~spin bands of V atoms crossing 
the Ep level, a result which we have also confirmecPH 
Our result is that inclusion of SOC in LSDA+U opens 
a band gap, signaling a half-metal-to-insulator transition 
as the SO coupling parameter is switched on. Since we 
argue that the uniform complex ferro order is the correct 
orbital order based on exchange constant calculations, 
we predict a half-metal-to-insulator transition to occur 
in single crystalline MnV204 as the temperature goes 
below Tg. In Fig.Qb) we present the corresponding band 
structure of the non-collinear magnetic configuration for 
this order, with the partial characters of V t<i g shown. We 
find that the Mn atoms carry no orbital moment, but the 
V atoms have an orbital moment m = 1.03. The spin 
moments are m s = 4.33 for Mn, and m s = 1.71 for V. 
Since the spin and orbital moments are antiparallel due 
to SOC coupling, the total moment for V is m w 0.7 in 
this phase. 

We note for comparison that Ref. H] also finds a uni- 
form (same) orbital occupied on all V sites, except that 
they relax the structure to obtain trigonal distortion of 
the V06 octahedra, and hence a 45° rotation of the sec- 
ond electron's t^g orbital on alternate V atoms of the 
afe-plane V chains. We do not have a trigonal distortion 
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in our setups, but we do find that the resulting exchange 
interactions correctly describe the spin-wave excitations 
of the low-T non-collinear phase (see below). This gives 
us a confidence that trigonal distortion has only a mini- 
mal influence on the magnetic exchange interactions, and 
that our LSDA+U+SO orbital order is equally good for 
the purpose of computing these interactions. For refer- 
ence, we label this order as 'SOC-FOO'. 



III. RESULTS FOR EXCHANGE 
INTERACTIONS 

Here we outline the spin-wave model, the ground state 
spin configuration, and present the results for our calcu- 
lated exchange constants J and single-site anisotropy pa- 
rameters D. Our obtained spin wave spectra of M11V2O4 
and comparisons with the neutron scattering experiments 



The subscripts on the J label the four inequivalent V 
sublattices, i, j, k, 1, and two inequivalent Mn sublattices, 
p, q. The Jmii-v is taken outside the summation because 
it has the same value for all pairs of Mn and V atoms. 
All the J couplings are between nearest-neighbor atoms 
of two different sublattices, and each pair is counted only 
once in the summation over all sites. We ignore the next- 
nearest-neighbor couplings because we found them to be 
much smaller in magnitude. 



B. Spin Configuration 

The low-T spin structure is non-collinear for V atoms, 
and collinear for Mn atoms with respect to the c-axis. 
The pyrochlore lattice on which the V atoms sit is 
geometrically frustrated for nearest-neighbor isotropic 
(J ab = J c ) AFM exchange. The frustrated pyrochlore 
interactions mean that there could be a macroscopic 
ground state degeneracy. But this frustration is par- 
tially relieved in the low-T phase by the presence of ad- 
ditional nearest -neighbor exchange interactions with Mn 
atoms, tetragonal distortion, and orbital ordering. The 
last one has the effect of making the V-V AFM exchange 
anisotropic: J ab ^ J c . It is well-known that the orbital 
or magnetic degeneracy can be lifted by the coupling of 
these degre es o f freedom with the lattice via the Jahn- 
Teller efFectpES. The ground-state spin configuration se- 
lected by the system in the low- T phase is non-collinear 
due to the combined effect of the frustration and coupling 



are also given. 



A. Spin Wave Model 

The parameters of the model are: [1] the exchange 
constants J derived from the LSDA+U(+SO) con- 
verged charge densities using linea r response theory and 
the magnetic force theorerri 22 * 23 *, and [2] the single- 
ion anisotropy parameters D calculated using an exact- 
diagonalization atomic multiplet procedure^. We input 
these parameters into the Heisenberg model Hamiltonian 
with anisotropy terms, minimize the classical energy to 
find the stable ground state configuration, and calculate 
the spin-wave excitation spectra. The model Hamilto- 
nian is: 

<jl) <kl> 

S q + Yl ^-D X -S K . (1) 

x=i,j,kj,p,q 
I 




FIG. 5. The T — non-collinear spin configuration of the V 
spins from Ref. |3| 

of the V spins to Mn spins, V t^g orbitals, and the lattice. 
In this structure, the V-atom spins develop components 
in the ab plane perpendicular to each other. The amount 
of canting away from the c-axis can be characterized by 
a single canting angle 8. Given the values of all J and D 
in Eq. [T] one can find the angle 8 as a function of J and 
D that will minimize the classical ground-state energy 



-Hspin — — ^ JijSi ■ Sj — JikSi ■ 5k — ^ JilSi ■ S\ — 22 ^jk<Sj • 5k — 
<ij> <ik> <il) <jk) 

(S P + S q ) ■ (3 + 5j + S k + §1) - X] J pq^p ' 



<(p,q)(i,j,k,l)) 



<pq> 
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of the configuration (derivation given in Ref. 



arccos 



3JMn-V<S , Mn 



(Dy ~ Dy - 2Jy_y - 2 Jy _y) Sy 



(2) 



The non-collinear spin configuratio n th at achieves this 
energy minimum is shown in Figure p^l 



C. Exchange Constants 

In Table [T] we present the J parameters that we calcu- 
late using LSDA+U method and magnetic force theorem. 
As the method computes the exchange constants in re- 
ciprocal space, we fourier transform them and show only 
nearest-neighbor exchange interactions between atoms of 
each sublattice. For the spinel structure , any of nearest- 
neighbor pairs always belongs to a different sublattice. 
The values of Jy_y f° r no orbital order and U = eV are 
the same for all V-V pairs; but when U — 5 eV, there is 
a tendency for anisotropy to develop: the in-plane Jy-v 



becomes unequal from the out-of-plane </y_y- This is 
one reason in support of our earlier argument that the 
anisotropy in the Jy_v, and the orbital ordering which 
causes it, are both interaction driven. When there is an 
orbital order, Jy_y is different along the ab V chains and 
between the chains (along the c-axis), as expected. This 
is true even for the case of the uniform orbital orders, 
because the exchange matrix elements of the Coulomb 
operator will be different within the ab plane and be- 
tween the planes, as can be seen from the shapes of the 
occupied orbitals in Fig. [l|a,b) . 



best match the experimentally known D values are given 
in Table HO for both the V 3d 2 and Mn 3d 5 shells. 

The easy axis for Mn is z(c)-axis, and for V it is either 
x or y. The easy axis always has a negative anisotropy 
parameter, which means the energy is lowered when the 
spin projection along the easy axis is maximized. For Mn, 
the spin projection along z tends to be maximized. How- 
ever, V also has a positive anisotropy parameter along 
the z-axis. So, V spin projection likes to be maximized 
along y, and minimized along Thus, the V spin mo- 
ment has a tendency to be in a non-collinear direction 
with respect to the z-axis. Our anisotropy computation 
is able to reproduce these signs as well as magnitudes for 
Mn (3d 5 ) and V (3d 2 ) shells. 

Looking at the anisotropy fit values we see that the 
value of Mn anisotropy reported in Ref. [15| are ob- 
tained for E xy = —0.016 eV, E eg = 1.0 eV, namely 
DJ, n = —0.1123 meV, which is similar to the literature 
value. For our fitted values of V anisotropy, we do not 
find such a close match, but there are several CF values 
which give the anisotropy of Ref.[T5lup to the correct sign 
and order of magnitude. For example, E xy — —0.024 eV, 
E eg = 0.4 eV give D v = 7.34 meV, and Dy' v = -4.056 
meV, which can be compared to D y = 2.79 meV, and 
D v ' v = —4.04 meV of Ref. [TS] One reason why our cal- 
culated Dy parameter differs by a large amount from 
the experimental value is because we have to tune the 
CF energy levels to simultaneously match two different 
single-ion anisotropies, and it was not possible to get 
them both to match the experimental D values of V. 



E. Spin— Wave Spectra 



D. Single— Ion Anisotropy 

The calculation of single-ion anisotropy requires first 
the total energies of the interacting atomic shell in a 
crystal-field environment along with the spin-orbit cou- 
pling. The method for its computation is described in 
Ref. [X5] and 1261 so here we merely present our results. 
The input parameters used in the total energy calcula- 
tion are: the SOC parameter 0.15 eV, and Slater inte- 
grals F = 5.0 eV, F 2 = 7.6 eV, F A = 4.7 eV. We then 
vary the direction of the magnetic moment by applying 
a small external magnetic field. The CF levels have the 
following values: the energies of yz and zx are both set to 
the reference value of 0.0 eV. The e g level is varied from 
0.2 eV to 1.0 eV in steps of 0.2 eV; and the xy level has 
the energies -0.024 eV, -0.016 eV, -0.008 eV, 0.0 eV. 
This gives a set of 20 different CF configurations. The 
E xy = E yz/xz represents cubic CF, and E xy ^ E yz / xz 
represents tetragonal CF. 

The total atomic shell energies thus obtained are fitted 
to a parabolic function of the polar angle 9 representing 
moment orientation, centered at 9 = in the case of z— 
axis anisotropy, and centered at 9 — ir/2 in the case of 
x/y-axis anisotropy. The results of the parabolic fit that 



We developed a code to compute the linear spin wave 
spectra for the non-collinear spin configuration. The pro- 
gram takes as an input our computed values of J and D. 
We first find the ground state which will in general be a 
non-collinear configuration with the spins pointing along 
the local quantization axis as given by 9 in Eq. [2j We sec- 
ond find the Heisenberg equations of motion, and numeri- 
cally diagonalize the resulting system of linear equations. 
The resulting spin-wave spectra are plotted in Fig. [6^ a) 
for the SOC uniform orbital order, along with the experi- 
mental spin waves. We find that the spin waves obtained 
from the J and D values of the SOC ferro-orbital order 
matches well with experiment. We also note that the 
lower two modes are due to the oscillations of Mn spins: 
The lower energy being the symmetric mode, and the 
higher energy the anti-symmetric mode^. The upper 
four modes are oscillations of the V spinal 

For comparison, we also show the spin waves for the 
other orbital orders do not match well with the experi- 
mental data. The model parameters for these orbital or- 
ders do give a reasonable spin canting angle when using 
Eq. |2j but the upper branches of the spin waves corre- 
sponding to the V oscillations are too high in energy and 
have a larger band-width in these plots (due to consid- 



meV 


No OO U = eV 


No OO [/ = 5 eV 


AFOO-I U = 5 eV 


AFOO-II [/ = 5 eV 


SOC-FOO U = 5 eV 


ExptP 


Ju 


-2.72 


0.136 


0.3264 


-0.04488 


0.1496 


- 


jab 

J V-V 


-20.4 


-21.76 


-14.96 


-19.04 


-11.56 


-9.89 


Jy-V 


-20.4 


-18.36 


-3.536 


-7.072 


-2.72 


-3.08 


JMn-V 


-10.2 


-2.992 


-5.44 


-5.44 


-4.76 


-2.82 


Jpq 


1.2 


2.167 


2.72 


2.72 


2.72 




Jpp 


-0.476 


0.204 


0.204 


0.272 


0.272 





TABLE I. Calculated exchange constants in meV for the collinear ferrimagnetic configurations and imposing various orbital 
orders. The last two columns list our theoretical Js for the spin-orbit coupled ferro-orbital order (SOC-FOO) and experimental 
J's from Ref. [15] respectively. The experimental J's match the SOC-FOO J's better than for the other theoretical J's. 






CF E xv eV 


CF E eg eV 


Theory meV 


ExptP me V 


Mn D z 


-0.016 


1.0 


-0.1123 


-0.1024 


V D x ' v 


-0.024 


0.4 


-4.056 


-4.04 


V D z 


-0.024 


0.4 


7.34 


2.79 



2.0 2.5 

->| 2tf(4, 0, q) 



FIG. 6. In all panels, the red lines are experimental and black 
lines are theoretical spin-waves: (a) Spin wave spectrum for 
I^i/amd spin-orbit coupled ferro-orbital order (SOC-FOO) 
along the high-symmetry lines of the Brillouin zone. We find 
an excellent match between our theoretical and previous ex- 
perimental data of Ref. 1151 (b) Spin-waves corresponding to 
the I4\/a symmetry AFOO I order. The upper four V os- 
cillation branches of the theoretical spin-waves are both too 
high in energy and have a larger dispersion compared to the 
experimental plot, (c) Same as in (b), but for the AFOO II 
order. Here the overestimate in the Jv-v is even greater than 
in (b). 



TABLE II. Table of calculated anisotropy constants for V 3+ 
3d 2 and Mn 2+ 3d 5 atomic shells. The energies of E xy and 
E Eg due to the crystal-fields are measured with respect to 
E vz , zx = eV. 



erable overestimate of V-V exchange) compared to the 
correct one in Fig.[6]Ja). Figure[6]Jb) shows the spin waves 
for the 'AFOO F order, which is composed of alternating 
yz and zx ab layers along the c axis. 

Figure |6jc) shows the spin waves for the 'AFOO IF, 
which is composed of real linear combinations of yz and 
zx orbitals with the relative sign between yz and zx al- 
ternating between ab layers along the c axis. As Table [I] 
shows, this order again gives a considerably greater V-V 
exchange compared to the experiment and therefore the 
upper branches are much higher in energy and have a 
greater dispersion relative to the experimental plot. 

We conclude that the excellent agreement between our 
theoretical and experimental spin-wave dispersions for 
all the six oscillation modes can be obtained only when 
using the exchange interactions derived from the SOC 
ferro-orbital order with IAi/amd symmetry, where the 
second ti g electron occupies a complex linear combina- 
tion of \yz) ± i\zx) uniformly on all V-sites. The other 
two orbital orders, 'AFOO I' and 'AFOO IF, do not give 
such a good match with experiment throughout the Bril- 
louin Zone, so these orders can be ruled out. We further 
note that the spin-orbit coupling plays an important role 
in the orbital physics of V-atoms in MnV204. This is 
also justified by the fact that the single-ion anisotropy 
is relatively high, as evidenced by the large gaps for the 
would-be acoustic modes at T. 



IV. CONCLUSION 

By theoretical computations of the interatomic ex- 
change constants using LSDA+U method, magnetic force 
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theorem and by imposing various orbital ordering scenar- 
ios we have shown that the orbital order on the V sites 
of MnV2 04 is similar to a complex linear combination 
of zx and yz on all V sites. Our calculated spin wave 
spectra for this order come closest to the experimental 
data. Further support in an evidence of the complex 
order is the strong single-ion anisotropy experienced by 
the spin moments on the V sites. We also predict, based 
on our U — 5 eV, orbital-ordered band-structures, that 
the low-T phase of MnV2 04 is a Mott-type insulator, 



and that a half-metal-to-insulator transition accompa- 
nies the simultaneous orbital ordering, structural distor- 
tion, and non-collinear moment transitions at T s = 53 
K. 
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